Strong electron-phonon coupling of the Fe breathing mode of LaOi-^F^FeAs 
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The electron-phonon coupling of LaOi-^F^FeAs is re-investigated on the basis of density func- 
tional theory in local density approximation. The implications of the nesting of the Fermi 
surfaces are carefully studied and found to lead to a non-standard electron-phonon coupling of the 
corresponding Fe in-plane breathing mode, which might make a strong contribution to the high su- 
perconducting transition temperature. The semi-metallic behavior of the undoped material is also 
further illuminated. 
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The recently discovered [l[ layered superconductor 
LaOi-^F^FcAs, 0.03 < x < 0.13, attracts at present con- 
siderable interest due to its comparably high critical tem- 
perature T c w 40K and possibly high upper critical field 
H c2 ~ 50Tesla, due to its similarities to the cuprates, and 
due to many unusual properties. One peculiarity of all 
density functional based electronic structure calculations 
12,11, II 0,001 is that, although there is strong Fe-3d As- 
Ap covalency resulting in a total width of the correspond- 
ing band complex of about 7eV, in a window of about 2.5 
eV around and below the Fermi level the electronic states 
are of nearly pure Fe-3d character (Fig. [T|) . Mazin et al. 
@ pointed out that there is a strong Fermi surface (FS) 
nesting in the undoped material with q = (7r, 7t, 0) nest- 
ing wave vector between the two hole cylinders around 
the tetragonal axis T — Z and the two electron cylin- 
ders around the axis M — A shifted by q against T — Z 
Fig\ [21. This nesting led those and several other authors 
1 la, 0] speculate on an important role of antiferro- 
magnetic correlations with magnetic in-plane superstruc- 
ture vector (it, tt), an apparently well justified scenario. 
Electron-phonon (e-p) coupling obtained on the basis of 
density functional perturbation theory [2, Q yielded only 
a weak coupling quite not enough to explain the high T c . 

However, as is well known 1(J 11 1 FS nesting also leads 
to phonon anomalies, and the corresponding e-p coupling 
may not be correctly described in lowest order perturba- 
tion theory which disregards the degeneracy. Moreover, 
since metallicity and superconductivity take place in the 
Fe layers of this quasi-two-dimensional tetragonal mate- 
rial only (Fig. [3J, an average of linearized e-p coupling 
over the whole structure might be misleading. In the 
present text we report results of investigation of these 
points which might be equally important compared to 
the magnetic correlations to understand the material. It 
is indeed found that e-p coupling may be quite strong. 

High precision density functional calculations were per- 
formed with the local density approximation functional 
of [IH using the full-potential local-orbital code [3 in 
the version FPL07-28 " 



FIG. 1: (color online) Total and component 
decomposed densities of states of LaOFeAs. 
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with its default orbital set- 
tings. The experimental lattice constants of undoped 
LaOFeAs, a = 4.03552A, c = 8.7393A, were used, 



FIG. 2: Bands close to the Fermi level (ep = 0) in the 
P4/nmm zone. 



and the two Wyckoff parameters were optimized to be 
z As = 0.365, z La = 0.1435. 

The phonon mode relevant in this respect is the in- 
plane breathing mode of the Fe square lattice (checker- 
board pattern where every other white square is shrunk, 
Fig. IU left part). A static displacement of Fe according 
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FIG. 3: (color online) The P4mm unit cell (Fe 
breathing mode) of LaOFeAs. 
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FIG. 4: Left:Displacement pattern of the Fe square 
lattice for the Fe breathing mode. The P4mm unit cell 
is shown (thick lines) and the P4/nmm unit cell 
(dashed lines, its center of inversion for the full 
structure is half-way between neighboring Fe sites). 
Right: Total energy per two formula units vs. Fe 
displacement of the Fe breathing mode. 



to this breathing mode reduces the P4/nmm symmetry of 
the lattice to P4mm by doubling the unit cell and remov- 
ing the inversion from the space group P4/mnm. (The 
center of inversion was half-way between two neighbor- 
ing Fe sites.) Doubling the unit cell of the lattice means 
folding down the unit cell of the reciprocal lattice, the 
Brillouin zone. The original point M : q = (ir, ir, 0) of the 
P4/nmm lattice (Fig. [2]) is folded onto r : q' = (0, 0, 0) of 
the P4mm lattice; the folded band structure from Fig.[2]is 
shown on Fig. [5] For symmetry reasons, at q = (n, tt, qs) 
the Fe-mode shown on Fig. |4] could only mix with a cor- 
responding oxygen breathing mode; both are, however, 
expected to interact only very weakly, and we neglect 
this mixing. The total energy change due to a static dis- 
placement (frozen phonon) is shown on Fig.[4j right part. 
Despite the FS nesting the mode potential is quite har- 
monic, the reason for this quadratic dependence on the 
Fe displacement uf c is given below. It corresponds to a 
frequency of wf c = 29.4meV in good agreement with the 
upper limit of the metal mode energies obtained in 
Note that there are two equivalent such modes of the Fe 
atoms as a rotation of the displacement vector by 90° in 
the left part of Fig. 0] corresponds just to an in-plane shift 
of the whole P4mm lattice by (1/2,0,0) which together 
with a reflection of the c-axis is a glide symmetry element 
of the original P4/nmm lattice. 

On Fig. [5] the bands corresponding to the undistorted 
lattice and to two static mode amplitudes, respectively, 
are shown in the downfoldcd P4mm zone. The common 
Fermi level is again £p = 0. Now, the nesting for q' = is 
seen as band crossings very close to the Fermi level (ob- 
serve the energy scale of the figure) on the lines T — M 
(and r — X) and Z — R (the latter for q' 3 = ir; the new 
point M of the folded zone corresponds to point X on 
Fig. [2]) From these bands as function of the displace- 
ment ufo of the Fe atoms, deformation potentials are 
derived. Care is, however, needed since deformation po- 
tentials arc touchy entities [l5j]. First of all, they have 
to be determined from a common Fermi level. In the 
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FIG. 5: (color online) Band deformations due to the Fe 
breathing mode in the P4mm zone. The upward 
dispersing band on the line T — M is two-fold 
degenerate. The inset on the right shows the zoomed 
detail from the line M — T. 

present case, moreover, reversion of up e is equivalent to 
an in-plane shift of the whole P4mm lattice by (1, 1,0) 
which is a lattice vector of the original P4/nmm lattice, 
and hence changes in band energies Ae;, Tl must be even in 
up c . Consequently, well away from band crossings there 
is no first order deformation potential, Ae^n ~ u Fo . At 
crossing of folded bands according to the symmetry of 
the lattice with uf c 7^ present (Fig. [5|) , gaps open with 
Ae ga p ~ |ufc|- In the vicinity of those crossings one has 

Ae kn « ±v/(«Afc) 2 + (IfcWFc) 2 T \vAk\, (1) 

where v is the average of the two band velocities at the 
crossing, Afe is the deviation from the crossing point, and 
IfcUfo is the deformation potential at the gap: 

Ae ga p « 2J F eMFc- (2) 

For the displacement corresponding to one phonon per 
unit cell with energy ajpc, ( u f c ) = {^Mfc^Fc)" 1 ^ 2 one 
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has Jpe^Fc) 1 / 2 ~ 20mcV~ 0.7wf c - By expanding the 
square root of ([1} for |t>Afc| = 20mcV, approximately 
the position of the FSs of the fluor doped material (at 
lower doping level compared to Fig. where 12.5 p.c. 
fluor were considered to limit the computational effort), 
an FS-averaged e-p coupling constant 



9 — 9kn.k' 



1 4e«) 

2 20meV 



10meV«/ 2 ,Fc«> (3) 



is obtained, where /2,Fe is the corresponding second or- 
der deformation potential. Of course, whenever a phonon 
mode connects two points on the FS, this gaping with 
lattice displacement happens. However, normally it hap- 
pens just for those two points without statistical weight. 
The peculiarity here is that, as a consequence of the nest- 
ing, for one mode it happens in the vicinity of the whole 
FS in a layer of thickness ~ w v in energy (cf. Fig. [S]). 

In the standard theory the e-p coupling strength A = 
Z(0) — 1 = i(9E/9e) e= o is obtained from the e-p self- 
energy 



„. , \ f duo f d 3 k 



; G n {e-Lo,k) (4) 



(n being the band index and v the phonon mode in- 
dex) by treating the phonon quantities as nearly k- 
independent and using the fc-integrated electron Green's 
function 



E 



d 3 k 



G n (e -u,k) 



N(0) / de c - 



~i(e — lo) — e c 
= -i7dV(O)sign(e-cj) + 0(l/e c ), (5) 

where -/V(O) is the electron density of states (DOS) at 
the Fermi level, and the high-energy cutoff e c is con- 
sidered negligible small (of the order of Migdal's pa- 
rameter). Then, the derivative of sign yields a delta 
function and the ^-integration is readily done yielding 
A = E„27V(0).g 2 M,. 

In the present case, which is not subject to Migdal's 
theorem, this approach is not allowed since the renormal- 
ized electron Green's functions have a low energy struc- 
ture at rjc « 20mcV, the distance of the gap edges from 
Fermi level, which causes additional terms in ([5]). More- 
over, Eq. ((3|) represents a two-phonon vertex (four- leg 
vertex), and the corresponding contribution to the elec- 
tron self-energy has two phonon propagators connecting 
the two vertices and an additional w-integration in the 
self-energy, and without the situation with the gaps one 
would again be left with A sa 'Yly N (Q)g 2 /w„, where only 
g is now the coupling strength j3| of the four-leg vertex, 
which in the present case is of the same order of magni- 
tude as the usual g of an ordinary three leg vertex. How- 
ever, due to the gap structure all frequency integrations 
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FIG. 6: (color online) Deformation of the total DOS due 
to the Fe breathing mode (states per two formula units). 



are entangled, and resonances of the type l/(w„ — rj c ) 
may appear in the course of self-consistency for G„ [16| . 
A quantitative theory for this case is complicated and 
not yet well developed. Non of the approximate assump- 
tions leading to (j4]) and (O can be maintained. So far 
one could only speculate on the magnitude of the effects. 

Now, the nesting is considered in more detail. A closer 
look on Fig. [5] reveals that a displacement ufo = 0.0645A, 
which does not much exceed the average displacement of 
one phonon per mode, already removes two of the four 
FSs of the in-plane bands (Fe^c^ and Fe-3<ij, z orbital 
character) by opening a gap of about 75meV. Only an 
electron FS of the shape of a cylindric shell around the 
tetragonal axis, quite isotropic perpendicular to it, sur- 
vives (and the hole pocked around point Z of the band of 
Fe-3d z 2 orbital character). Just for illustration the effect 
of a (not quite realistic) two times larger displacement 
is also plotted on Fig. [5] for which the material becomes 
a semi-metal with only two tiny FS pockets, an electron 
ring and a hole pocket, both around Z. Even slightly 
larger displacements would make it a true semiconduc- 
tor. This is reflected in the behavior of the total DOS 
shown on Fig [6] (where the normalization is for the unit 
cell of the undistorted P4/nmm structure). One impor- 
tant point is that the nesting is only approximate and, 
more importantly, that at a given energy the bands which 
were originally electron or hole bands move in the same 
direction. Since the band crossings are not exactly at 
Fermi level and in e-p coupling theory the Fermi level 
has to be kept fixed [la ], the considered c-p coupling is 
transmitted to all bands and is also transmitted to the 
band of Fe-3d z 2 character. These findings (together with 
the similar magnetic peculiarities Q) are capable of ex- 
plaining many of the unusual electronic properties of the 
undoped material. 

Finally, the effect of fluor doping for oxygen is studied 
by replacing one of eight oxygen atoms with a fluor atom 
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FIG. 7: Bands of La0 .875Fo.i25FeAs in the Pmm2 
zone. 
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LaOi-^F^FcAs is indeed very strong due to the Fe in- 
plane breathing mode, and this cannot be treated in stan- 
dard e-p theory of superconductivity. With respect to 
this breathing mode, the material rather compares with 
MgB2, although the details are quite different, and the 
magnetic couplings already discussed in the literature 
must also be taken seriously into account. In that re- 
spect, the material may have much in common with the 
cuprates. A principally similar but in detail quite differ- 
ent case is a breathing mode of the copper square lattice 
in connection with the approximate (vr, 7r)-nesting [l6j |. 
It is intriguing hence, whether LaOFeAs could provide 
some bridge between the other two cases, MgB2 and the 
cuprates. The undoped material appears to be a bad 
metal with a great complexity of very strong couplings, 
e-p and magnetic. Doping (potentially both electron or 
hole) moves LaOFeAs away from this messy coupling 
situation into a still strong coupled but good metallic 
regime. It is to be expected from former Q and the 
present analysis that LaOFeAs will be more symmetric 
with respect to both electron and hole doping than the 
cuprates. 

Helpful discussion with Klaus Kocpcrnik, Stcfan- 
Ludwig Drcchsler, B. Biichner and in particular with M. 
Eschrig are gratefully acknowledged. I.I. Mazin pointed 
out to me a substantial error in an earlier version of the 
paper. 



FIG. 8: Total DOS of ordered LaO . 8 75F .i25FeAs 
(states per two formula units). 



in an ordered way. This reduces the symmetry of the 
lattice further to the orthorhombic Pmm2 group with a 
four times larger unit cell, square in plane, compared to 
the P4/nmm structure. The corresponding bands (again 
folded, into the Pmm2 zone, which however does not re- 
sults in new shifts of the FSs which are already close 
enough to the central line T — Z) and total DOS are 
shown on Figs. [7] and [H As can be seen by comparing 
Figs. [5] and [6] with Figs. [7] and [8j one fluor atom do- 
nates just one electron without changing the Fc derived 
^2g-bands (d xz , d yz ) in the vicinity of the Fermi level. 
Therefore, also no disorder effect of fluor doping is ex- 
pected in the relevant electronic structure. Moreover, the 
Fe-3d z 2 derived bands, which shift a little bit due to the 
change of layer charges, are now completely filled (one 
of them producing the DOS peak at about -150meV), 
and the corresponding hole pocket around Z has gone, 
leaving only the four slightly deformed cylindrical FSs, 
two larger electron cylinders around the T — Z axis and 
two smaller hole cylinders around the M — A axis of the 
original P4/NMM structure (cf. Fig. [2]). 

In summary, we have shown that the e-p coupling in 
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